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performed  by  Major  George  D.  Ballentine,  Air  Force  Weapons  Laboratory,  Kirtland 
Air  Force  Base,  New  Mexico. 
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SECTION  I 
INTRODUCTION 

The  Air  Force  must  periodically  evaluate  airfield  pavements  to  establish 
their  operational  capabilities  and  program  any  required  maintenance.  A  new  method 
of  nondestructively  evaluating  the  condition  and  load  carrying  capacity  of  air¬ 
field  pavements  is  currently  under  development,  but  one  shortcoming  of  that 
method  is  that  layer  thickness  information  is  required  to  evaluate  the  pavement. 

At  the  present  time  there  is  no  method  available  to  nondestructively  determine 
the  thickness  of  the  various  layers  of  a  pavement  system,  thus  necessitating  the 
use  of  "as-built"  drawings  or  digging  a  test  pit. 

This  study  was  undertaken  to  determine  the  feasibility  or  using  active  micro- 
wave  techniques  to  differentiate  between  and  accurately  measure  the  thickness  of 
various  layers  in  a  pavement  system.  Some  previous  work  has  been  done  to  measure 
the  thickness  of  Portland  Cement  Concrete  (PCC)  and  Asphaltic  Cement  Concrete 
(ACC)  materials  with  relatively  good  results,  but  very  little  is  known  about  the 
applicability  of  this  technique  to  construction  materials  below  the  PCC  or  ACC 
surface  course. 

When  using  active  microwave  techniques  to  characterize  these  natural  con¬ 
struction  materials  three  factors  are  of  prime  importance:  soil  type,  moisture 
content,  and  the  frequency  of  the  microwave  signal.  The  first  step  in  demon¬ 
strating  the  ability  to  measure  the  thickness  of  these  material  layers  is  to  show 
the  ability  of  the  microwave  signal  to  penetrate  the  material  under  varied 
moisture  contents.  The  measure  of  this  ability  to  penetrate  used  in  this  report 
is  the  skin  depth  or  e-folding  distance,  which  is  the  depth  at  which  the  wave  has 
been  attenuated  to  1/e  of  its  original  value.  (The  teorietical  model  which  is  the 
basis  of  this  work  is  included  as  an  appendix.)  If  the  skin  depth  is  sufficient 
to  indicate  that  the  wave  can  penetrate  the  full  depth  of  the  pavement  system 
(which  varies  from  61  to  122  cm),  and  that  for  similar  materials  under  similar 
conditions  the  measured  skin  depth  is  the  same,  the  feasibility  of  using  this 
technique  to  measure  pavement  thickness  would  be  established. 
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SECTION  II 

EXPERIMENTAL  PROCEDURE 

The  testing  in  this  study  involved  measuring  the  skin  depth  in  two  different 
soils  (gravel  and  clay),  at  three  different  saturation  levels  in  each  (10,  50, 
and  90  percent),  at  microwave  frequencies  in  the  range  of  0.5  to  4.5  GHZ. 

The  gravel  was  blended  by  mixing  concrete  aggregates  in  calculated  propor¬ 
tions,  so  that  the  grain-size  distribution  curve  as  shown  in  figure  1  was  ob¬ 
tained.  Soil  so  prepared  classifies  as  a  poorly  graded  gravel  (GP)  in  the 
Unified  Soil  Classification  system  (see  ref.  4).  Compaction  characteristics 
were  obtained  by  the  modified  Proctor  technique  and  are  given  in  figure  2.  From 
figure  2,  a  moisture  content  was  selected  for  making  soil  for  the  test  bed  in 
order  to  achieve  a  compacted  degree  of  saturation  of  approximately  10  percent. 

The  actual  degree  of  saturation  of  the  test  soil  (based  on  volume  calculations) 
could  not  be  accurately  predicted  until  the  bed  had  been  prepared. 

Two  boxes  of  identical  dimensions  (1.33  m  x  1.22  m  x  20.3  cm)  were  built  for 
making  two  beds  of  soil  simultaneously.  Slight  variation  resulted  in  the  density 
of  the  soil  bed  prepared.  The  boxes  were  lined  with  polyethylene  plastic  sheets 
to  prevent  leaking  of  soil  moisture. 

Mixing  of  the  dry  gravel  with  the  water  was  carried  out  in  a  concrete  mixer, 
and  compaction  of  the  bed  was  accomplished  by  means  of  a  gasoline-powered  vibra¬ 
tory  machine.  After  compaction,  the  box  was  weighed  by  supporting  it  on  a  knife 
edge  on  one  side  and  a  proving  ring  on  the  opposite  side.  Thus,  with  the  wet 
density,  moisture  content  and  grain  specific  gravity  known  for  the  soil  bed,  the 
void  ratio,  dry  density,  and  the  degree  of  saturation  could  be  calculated. 

The  above  techniques  were  employed  in  making  the  soil  beds  at  10  percent 
degree  of  saturation.  After  measurements  had  been  made  on  them,  the  degree  of 
saturation  of  the  soil  was  brought  up  to  50  percent  by  introducing  water  from 
the  surface  in  an  amount  calculated  from  the  dry  density  of  the  prepared  soil. 

The  box  was  again  weighed  to  provide  data  for  accurate  determination  of  new  wet 
density,  and  therefore,  degree  of  saturation.  Thus,  by  this  method,  compacted 
dry  density  of  the  soil  stayed  constant  from  one  saturation  state  of  the  test 
specimen  to  another.  The  same  technique  was  also  used  to  bring  saturation  to 
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Figure  2.  Modified  Proctor  Compaction  of  Gravel 
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90  percent,  although  the  preserve  of  trapped  air  in  the  soil  presented  some 
difficulties  in  achieving  that  high  degree  of  saturation. 

Data  on  the  two  soil  beds  at  each  of  the  three  saturations  are  given  in 
table  1.  Target  values  for  moisture  content  were  10,  50,  and  90  percent.  Table 
1  gives  the  act"al  values  for  each  sample. 

Some  difficulties  were  encountered  in  locating  a  soil  with  a  CL  classifica¬ 
tion  because  of  the  large  quantities  required.  Finally,  a  soil  was  selected 
which  had  to  be  dried,  crushed,  and  sieved.,  retaining  the  minus  number  50  sieve 
portion.  Soil  used  in  the  experiments  has  the  grain-size  distribution  curve  as 
shown  in  figure  3,  which  also  gives  its  Atterberg  limits.  Thus,  it  classifies 
as  a  CL  in  the  United  Soil  Classification  System.  Compaction  characteristics 
are  shown  in  figure  4. 

Dry  soil  was  mixed  with  an  appropriate  amount  of  water  in  a  soil  mixer  in 
batches  of  110  kg.  After  placing  the  soil  in  the  box,  it  was  necessary  to  compact 
it  with  a  hand  compactor  to  achieve  the  desired  density  and,  therefore,  the 
required  degree  of  saturation.  After  measurements  were  taken  of  this  soil  bed, 
the  soil  was  remixed  with  additional  water  to  bring  it  up  to  the  moisture  content 
necessary  to  achieve  50  percent  saturation.  This  time,  however,  a  gasoline- 
powered,  rammer-type  compaction  machine  was  used  to  attain  the  higher  density. 

For  the  90  percent  saturation,  the  soil  bed  with  50  percent  saturation  was 
percolated  with  water  through  four  plumbing  fittings  on  the  bottom  of  the  box. 

This  process  of  increasing  the  saturation,  aimed  at  maintaining  a  reasonable 
compaction  of  the  soil  (that  of  the  50  percent  saturated  sample),  and  avoiding 
the  need  to  compact  a  clay  soil  at  90  percent  saturation,  took  about  four  days 
to  accomplish  and  worked  satisfactorily.  Data  for  the  clay  beds  are  shown  in 
table  1 . 
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Table  1 


PROPERTIES  OF  SOIL  BEDS 


Soil  Bed 

Moisture 

Content 

(percent) 

Dry 

Density 

(g/cmJ) 

Degree  of 
Saturation 
(percent) 

Gravel  (Bed  1 ) 

2.0 

1.802 

11.4 

8.6 

1.802 

48.9 

15.1 

1.802 

85.9 

Gravel  (Bed  2) 

2.0 

*3- 
1 — 
CO 
• 

12.8 

7.4 

1.874 

47.3 

12.7 

1.874 

81.2 

Clay 

3.5 

1.352 

9.6 

11.0 

1.528 

39.3 

23.2 

1.528 

82.8 

10 
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Figure  2.  Modified  Proctor  Compaction  of  Gravel 


DENSITY,  g/cm 


Figure  4.  Modified  Proctor  Compaction  of  Clay 
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SECTION  III 
EXPERIMENTAL  RESULTS 

Preliminary  experiments  showed  that  by  sealing  the  samples  in  polyethylene 
the  moisture  was  retained  in  the  sample  in  a  fairly  uniform  state.  A  measure  of 
this  uniformity  was  obtained  by  comparing  surface  dielectric  constant  data  with 
the  dielectric  constant  for  the  layer. 

The  measurements  were  made  with  a  commercial  microwave  network  analyzer  with 
an  accuracy  of  +0.5  db.  This  particular  instrument  would  be  expected  to  produce 
the  greatest  skin  depth  uncertainty  at  the  lower  frequencies  and  saturation 
values.  Table  2  contains  information  on  the  sample  thickness  used  for  testing 
the  various  samples  at  the  three  different  frequencies.  The  thin  samples  at  high 
saturation  and  high  frequency  were  necessitated  because  of  the  high  attenuation 
of  the  signal  in  these  materials. 


Table  2 

SAMPLE  THICKNESS 


Soil 

Percent 

Saturation 

Frequency 

(GHz) 

Thickness 

(cm) 

Gravel 

10 

0.5 

20.2 

2.5 

20.2 

4.5 

20.2 

50 

0.5 

20.2 

2.5 

20.2 

4.5 

7.6 

90 

0.5 

20.2 

2.5 

20.2 

4.5 

9.7 

Clay 

10 

0.5 

20.2 

2.5 

20.2 

4.5 

10.8 

50 

0.5 

8.9 

2.5 

8.9 

4.5 

8.9 

90 

0.5 

2.5 

2.5 

2.5 

4.5 

2.5 
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The  skin  depths  calculated  from  the  measured  data  are  shown  in  table  3.  A 
statistical  analysis  of  this  data  was  performed  in  an  attempt  to  quantify  some 
preliminary  observations.  Because  the  experiment  performed  was  a  modified  23 
factorial  design,  it  produced  more  data  than  actually  needed  for  a  23  statistical 
analysis.  To  make  the  calculations  as  simple  as  possible  and  still  retain  data 
representative  of  actual  field  conditions,  analysis  was  performed  on  two  sets  of 


iata. 

The  first  set  of  data 

follows: 

Type 

Percent 

Frequency 

Skin  Depths 

Soil 

Saturation 

(GHz) 

Replication  1 

Replication  2 

Clay 

10 

0.5 

35.8 

33.1 

2.5 

19.3 

13.8 

50 

0.5 

10.1 

12.9 

2.5 

8.4 

5.6 

Gravel 

10 

0.5 

507.2 

980.0 

2.5 

77.9 

73.2 

50 

0.5 

304.7 

196.1 

2.5 

19.8 

13.3 

Treating  the  soil  type  as  factor  A,  the  degree  of  saturation  as  factor  B, 
and  the  frequency  as  factor  C,  an  analysis  of  variance  was  performed  on  these 
data,  with  the  following  results: 


Analysis  of  Variance 

D.F 

S.S. 

M.S. 

F 

Replications 

1 

7475.26 

7475.26 

... 

Treatments 

7 

912827.06 

130403.87 

8.28* 

Error 

7 

110202.24 

15743.18 

— 

*  Significant  at  the  1 

percent 

level 

The  individual  factor 

effects 

were  further  analyzed 

as  follow: 

Factor 

S.S. 

F 

A 

258495.98 

16.42* 

B 

85395.45 

5.42* 

AB 

67717.05 

4.30 

C 

213698.18 

13.57** 

AC 

193050.39 

12.26* 

BC 

50254.43 

3.19 

ABC 

44215.58 

2.81 

*  Significant  at  the  1  percent  level 

**  Barely  misses  being  significant  at  the  5  percent  level 
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Table  3 

MEASURED  SKIN  DEPTHS 


Type 

Soil 

Percent 

Saturation 

Frequency 

(GHz) 

Skin  Depths  (cm) 

Replication  1 

Replication  2 

Clay 

10 

0.5 

35.8 

33.1 

2.5 

19.3 

12.8 

4.5 

10.0 

11.2 

50 

0.5 

10.1 

12.9 

2.5 

8.4 

5.6 

4.5 

6.3 

8.5 

90 

0.5 

5.1 

6.2 

2.5 

4.1 

3.4 

4.5 

2.0 

2.0 

Gravel 

10 

0.5 

507.2 

980.0 

2.5 

77.9 

73.2 

4.5 

25.0 

17.3 

50 

0.5 

304.7 

196.1 

2.5 

19.8 

13.3 

4.5 

6.0 

4.6 

90 

0.5 

57.2 

44.1 

2.5 

5.6 

9.0 

4.5 

2.5 

2.5 

As  a  further  check  on  the  validity  of 
was  selected  for  a  similar  23  statistical 

used  in  this  set  of  calculations: 

Type  Frequency 

Soil  Saturation  (GHz) 

the  data  analysis,  a  second  set  of  data 
analysis.  The  data  shown  below  were 

Skin  Depths 

Replication  1  Replication  2 

Clay 

50 

0.5 

10.1 

12.9 

2.5 

8.4 

5.6 

90 

0.5 

5.1 

6.2 

2.5 

4.1 

3.4 

Gravel 

50 

0.5 

304.7 

196.1 

2.5 

19.8 

13.3 

90 

0.5 

57.2 

44.1 

2.5 

5.6 

9.0 
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Again  treating  the  soil  type  as  factor  A,  tiie  degree  of  saturation  as  factor 
B,  and  the  frequency  as  factor  C,  an  analysis  of  variance  was  performed  on  these 
data,  with  the  following  results: 


Analysis  of  Variance 

D.F. 

S.S 

M.S. 

F. 

Replications 

1 

967.21 

967.21 

— 

Treatments 

7 

100523.06 

14360.44 

19.90* 

Error 

7 

5051.17 

721.60 

— 

*  Significant  at  the  1 

percent  level 

The  individual  factor 

effects 

were  analyzed  as 

follows: 

Factor 

S.S. 

F. 

A 

22052.24 

30.56* 

B 

11891.92 

16.48* 

AB 

9990.00 

13.84* 

C 

20107.24 

27.86* 

AC 

18333.16 

25.41* 

BC 

9321.92 

12.92* 

ABC 

8826.60 

12.23** 

*  Significant  at  the  1  percent  level 
**  Significant  at  the  5  percent  level 

In  each  test  it  is  quite  clear  that  the  soil  type  is  significant  at  the 
1  percent  level,  thus  indicating  that  the  skin  depth  is  highly  dependent  on  the 
soil  type.  However,  it  is  also  apparent  that  the  moisture  content  is  highly 
significant,  thus  indicating  that  the  moisture  content  of  the  various  layers  of 
the  pavement  system  would  have  to  be  known  to  determine  their  thicknesses  by  this 
method.  The  frequency  of  the  microwave  signal  is  also  a  significant  factor  in 
the  test,  and  much  greater  skin  depths  are  measured  at  the  lower  frequencies. 

The  interaction  effects  are  also  quite  significant,  especially  at  the  higher 
moisture  contents. 

In  addition  to  the  above  statistical  analysis,  it  is  quite  apparent  from 
looking  at  the  data  in  table  3  that  there  is  a  huge  variation  from  replication  1 
to  replication  2  in  some  of  the  gravel  samples.  This  amount  of  variation  makes 
the  possibility  of  accurate  determination  of  layer  thickness  in  gravel  very  remote 
using  this  technique. 
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SECTION  IV 


CONCLUSIONS 

The  objective  of  this  effort  was  to  investigate  relationships  that  exist 
between  microwave  penetration,  soil  type,  moisture  content,  and  microwave  fre¬ 
quency,  and  thus  determine  if  it  is  feasible  to  use  microwave  techniques  to 
accurately  measure  the  thickness  of  various  soil  layers  in  a  pavement  system. 

The  analysis  of  data  in  section  III  clearly  indicates  that  the  moisture 
content  plays  a  significant  role  in  the  attenuation  of  microwave  signals.  There¬ 
fore,  it  is  necessary  to  know  the  moisture  content  of  each  layer  and  its  effect 
on  attenuation  within  that  soil  before  the  thickness  of  that  layer  can  be 
determined.  There  is  presently  no  way  to  nondestructively  determine  the  moisture 
content  of  the  individual  layers,  so  there  appears  to  be  no  way  to  isolate  the 
moisture  factor  and  make  the  microwave  system  work. 

There  are  several  aspects  of  equipment  selection  and  usage  which  could  be 
considered  to  improve  the  accuracy  of  the  data  collected.  Mismatch  between 
antennas  and  the  measuring  system  could  produce  a  systematic  error  which  could 
be  minimized  by  using  VSWR  antennas.  The  accuracy  of  the  microwave  network 
analyzer  could  be  improved,  and  additional  work  could  be  done  to  confirm  the 
assumption  of  single  ray  theory  used  in  these  measurements.  It  is  unlikely, 
however,  that  any  of  these  steps  could  overcome  significance  of  the  moisture 
content  on  the  measurements. 

Because  of  the  statistical  significance  of  the  moisture  content  and  the 
interactive  effects,  and  the  complexity  of  isolating  these  effects,  it  is  recom¬ 
mended  that  no  further  research  be  performed  in  an  attempt  to  use  microwave 
techniques  to  measure  pavement  layer  thicknesses. 
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APPENDIX 

THEORETICAL  MODEL 


The  following  equations  outline  the  properties  of  a  complex  medium  (see 
ref.  1).  In  general,  wave  propagating  in  a  medium  undergo  attenucition  and  phase 
shift  as  they  propagate.  The  complex  propagation  constant  for  plane  waves  is 

y  =  jw  /me  |^l-j  (1) 

and  can  be  written  as 

Y  =  a(w)  +  jB(w)  (2) 


Solving  for  a(w)  and  B(w)  gives 


a(cj) 


8(w) 


which  for  small  losses  become 


a (w)  ~ 


8  w  ~  oj  /  ye 

The  intrinsic  impedance  of  a  medium  is  given  by 


(3) 


(4) 


(5) 

(6) 
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The  model  used  consisted  of  one  layer  of  finite  thickness,  of  the  material 
under  test  sandwiched  between  air  and  an  aluminum  sheet  (see  figure  5).  The 
medium  is  assumed  to  be  isotropic  and  homogeneous  with  plane  boundaries.  Elec¬ 
trical  properties  o*  the  layer  are  described  by  permeability  (y),  permittivity 
(e),  and  conductivity  (o).  The  layer  is  assumed  nonmagnetic  (y  =  y0,  where 
Uo  *  4  (10) 7  hy/m).  Additionally,  e  is  assumed  to  be  Independent  of  frequency. 
Ray  paths  for  electromagnetic  energy  are  also  assumed. 

Due  to  the  directivity  of  the  ideal  horns  assumed  in  the  model,  only  reflec¬ 
ted  waves  are  received.  In  practice,  considerable  design  placement  and  polariza¬ 
tion  strategy  are  required  to  reduce  the  direct  coupling  to  an  acceptable  level. 
The  horns  used  were  not  ideal,  however  a  single  ray  normal  Incidence  model  was 
used  here. 


The  input  reflection  coefficient  in  medium  1  for  a  single  layer  can  be 
written  as  (see  ref.  2) 


-  e‘^2^2 

1  - 


(8) 
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where  r ,  is  the  reflection  coefficient  of  a  semi-infinite  medium  and  is  approxl 
mately  given  by 


ri 


(9) 


and  y2  by  (2),  (3),  and  (4).  The  magnitude  of  r.n  will  go  through  successive 
maxima  and  minima  as  a  function  of  frequency,  layer  thickness,  dielectric  con¬ 
stant  and  loss.  Minima  (or  maxima)  will  occur  approximately  every  2tt  radians. 
Thus,  when  gj  is  varied,  successive  minima  occur  when 


2UM  =  2tt 

from  which 


(10) 

(ii) 


where  successive  minima  occur  at  frequencies  fj  and  f2. 

It  has  been  shown  (see  ref.  2)  that  an  approximate  expression  for  the  maxima 
of  rin  can  be  obtained  when  cos  262£2  =  1.  Equation  (8)  then  becomes 


r  -2a2Jl2 

_  rt  -  e  22 
m  ~  1  „  -2a,  £ 


1  -  rie 


2  2 


Solving  for  the  attenuator  factor 


-2a2Jt2 


(12) 


(13) 


The  reflection  r,  as  given  in  (9)  will  always  be  negative,  thus  the  coeffi¬ 
cient  rm  as  given  in  (12)  is  also  negative.  Because  |rj  >  J rx  | ,  logarithms  of 
numbers  are  involved. 

Figure  6  shows  how  the  reflection  coefficients  in  this  model  are  related. 

The  jrj  and  |rin|  are  deduced  from  measurements.  The  |rj  is  determined  by 
using  frequencies  f!  and  f2  in  (11).  The  | r ^ n |  is  determined  by  return  loss 
measurement  techniques  (see  ref.  3).  The  sample  thickness  is  held  constant  thus 
the  attenuation  of  the  sample  can  be  determined  from  (13). 
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Reflection  from  reference  plate  on  top  of  sample 
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Figure  6.  Reflection  Coefficient  Relationship 
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Skin  depth  Is  defined  as  the  depth  at  which  the  wave  has  been  attenuated  to 
1/e  of  Its  original  value.  Thus,  skin  depth,  6,  can  be  determined  from  the 
relation 


a?6  =  1 


(14) 


ABBREVIATIONS  AND  SYMBOLS 


Signal  level  with  respect  to  1  mW,  decibels 

Incident  electric  intensity,  v/m 

Reflected  electric  intensity,  v/m 

Grain  specific  gravity 

Degree  of  saturation 

Velocity  of  light  in  free  space,  m/sec 

Natural  logarithm  base 

Frequency,  Hz 

1,  2  for  media  1  and  2  respectively 
Sample  depth,  m 
Time,  sec. 

Reflection  coefficient  of  infinitely  thick  sample 

Reflection  coefficient  of  finite  thickness  sample 

Maximum  reflection  coefficient  of  finite  thickness  sampl 

Attenuation  constant,  nepers/m 

Phase  constant,  radians/m 

Propagation  constant 

Dry  density,  g/cm 

Wet  density,  g/cm 

Skin  depth,  m 

Permittivity,  farad/m 

Permittivity  relative  to  free  space,  farad/m 

Permeabil ity,  henry/m 

Conductivity,  mho/m 

Angular  velocity,  radians/sec 

Intrinsic  impedance,  ohm 
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